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Abstract 

We review the experimental constraints on additional macroscopic Yukawa 

forces for interaction ranges below 1 cm, and summarize several theoretical 

predictions of new forces in this region. An experiment using 1 kHz mechanical 

oscillators as test masses should be sensitive to much of the parameter space 

covered by the predictions. 
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I. INTRODUCTION 



Experimental searches for a deviation from Newton's law of universal gravitation have 
received a great deal of attention over the past three decades. Motivation has come from 
reports of experimental anomalies, and from new theoretical predictions. Searches have 
involved a great variety of experiments on the planetary, lunar, geological, and laboratory 
scales. A recent review by Fischbach and Talmadge []I| shows that new forces with a strength 
weaker than or comparable to gravity have been excluded over distances ranging between 
1 and 10^^ cm. Experiments have only marginally explored the distance range under 1 
cm, however, and there is little knowledge of gravity itself in this range 0. Furthermore, 
the sub-centimeter region is of increasing experimental interest, given a number of recent 
predictions of new forces from modern theories which attempt to unify gravity with the 
strong and electroweak interactions. 

In this paper we present a review of the existing experimental limits and some of the 
recent theoretical predictions. We then describe an experiment which we expect to be 
sensitive to gravity at distances as short as 50 /xm, and to much of the parameter space for 
the new predictions. 

II. CURRENT LIMITS 

Experimental results from searches for new long-range forces of gravitational strength 
have been parameterized in several ways The Yukawa interaction is most commonly used, 
and we adopt it here. The potential due to gravity and an additional Yukawa interaction 
between two macroscopic objects may be written 

V{r) = - fdn /t/r2^^i^^^i^^^[l + aexp(-ri2/A)], (1) 
J J ri2 

where G is the Newtonian gravitational constant, is the distance between ri and r2, pi{ri) 
and P2{f^2) are mass densities, a is the strength of the Yukawa force relative to gravity, and 
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A is the range. If, as is likely, the Yukawa forces couple to quantities other than mass, then 
a will have a slight composition dependence which we do not consider. 

The current experimental limits on a, for Yukawa ranges between 1 /xm and 1 cm, are 
illustrated in Fig. |I], together with the limit anticipated from the experiment outlined in 
Sec. |V| below. Also indicated are specific predictions for new long-range forces, described 
in Sec. |T|. The region above and to the right of the bold curves is excluded by experi- 
ments falling into two general categories: classical gravity measurements and Casimir force 
measurements. 



A. Classical Gravity Experiments 

For values of A greater than lOO/im, the most stringent limits on a can be obtained 
from torsion balance experiments designed to test the gravitational inverse-square law. Well 
before the interest in gravitational measurements brought on by suggestions of a "fifth force," 
the report of a possible deviation in a measurement by Long motivated several similar 
experimental programs (see, for example, Refs. |0-§|) which have apparently ruled out any 
evidence for the anomaly [^. 



One such experiment is the Irvine 2-5 cm null measurement [|I0|, which yields the best 
limits in the range near 1 cm. This experiment used a torsion balance to measure the 
force on a test mass suspended inside a hollow cylinder, a non-zero result for which would 
have been evidence for an anomalous force. After an exhaustive consideration of possible 
systematic effects, and careful analysis of inhomogeneity in the cylinder, the authors derived 
the limit curve reproduced in Fig. |I], corresponding to a force sensitivity of about 4 x 10^^° 
dyne. Only the curve corresponding to the (more sensitive) limiting case of a hypothesized 
attractive Yukawa force [a > 0) is shown. 

Below the range of a few millimeters, the most stringent limits may be obtained from 



the torsion balance experiment of Mitrofanov and Ponomareva []TT|. This experiment used 



a dynamic-resonant technique to distinguish the gravitational interaction of the test masses 
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from backgrounds. A test mass, initially spaced 3.8 mm perpendicular to one end of the 
balance, was moved synchronously with the oscillations of the balance. The authors were 
able to achieve a sensitivity of about 3 x 10^^° dyne, for integration times of about 10'^ s. 
The ratio of the forces measured at 3.8 mm and 6.5 mm was used to obtain a curve similar 
to that shown in Fig. |I[ 



We have re-analyzed the data of this experiment, using Eq. 7 of Ref. |TT|, to extend the 
limit to shorter distances than those considered by the authors. The resulting curve is shown 
in Fig. ||. The la errors on the test mass separations, apparently not taken into account in 
the original report, are added in quadrature to the la errors on the measured force ratio, 



so that the limit we derive is slightly less sensitive than that reported in Ref. |TT]. We have 
not illustrated the cases for the separate hypotheses of attractive and repulsive forces, which 
have little effect on the results. 

Other inverse-square law and Cavendish experiments have been performed with inter- 



mass separations on the order of 1 cm, including that of the Tokyo group ||T2[. None, 
however, reach stronger limits in the range below 1 cm. The Rot- Wash experiment, which 
uses a 2620 kg depleted uranium attractor rotated around a torsion balance with lead and 



copper test masses, yields a limit nearly competitive with the Irvine result at 1 cm |TB|. This 
assumes a Yukawa force which couples to mass; the use of uranium makes this experiment 
especially sensitive to composition-dependent forces. Assuming a coupling to N — Z, the 
Rot- Wash experiment has set a limit 10 times more stringent than the Irvine experiment at 
1 cm. 

B. Casimir Force Experiments 

In 1948, H. G. B. Casimir showed that a weak attractive force should be present be- 
tween a set of flat, parallel conducting plates, due to the zero-point fluctuations of the 
electromagnetic field |l^. This force is given by 



Fc = —^A, (2) 
' 240 r4 ' ^ ^ 
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where r is the separation and A is the area. A similar effect, called the retarded Van der 
Waals force, occurs with dielectric plates and is given by an expression which differs from 
Eq. ^ only by a numerical factor of order unity [|T^ . 

For a given value of A, the most stringent limits on a can be obtained from Casimir 
and Van der Waals force data taken with mass separations of order A, given the strong 
dependence of Eq. |^ on r. These experiments usually investigate separation distances well 
below 1 mm. In practice, the best limits to date are derived from measurements of the force 
between a flat plate and a spherical shell, which varies as r~^. The experiment by Derjaguin, 
Abrikosova, and Lifshitz [0 done in 1956 provided the most sensitive limits in the range 



less than 100 /im for 40 years, until the recent measurement by Lamoreaux ||17|| . 

This experiment measured the attraction between a fiat, metal-plated quartz disk and a 
metal-plated spherical lens separated by distances ranging from about 0.5 to 10 fim. Using 
data from this experiment, we have derived the limits presented in Fig. |I[ For several values 
of A in the range 1-100 /im, a least-squares fit is made with the data in the measured force 
plot (Fig. 3, top, of Ref. ^7^, to the expression 

Fl^ = Fj{ai + ao) + ^ + h + Fy [a. A, p^, dj, {at + ao)]. (3) 

Oj -\- (Xq 

This is based on Eq. 7 of Ref. |l^, where the first three terms are the Casimir, electromagnetic, 
and constant terms used by Lamoreaux to fit the data (with the same fit parameters ao, /3, 
and h). We have added a Yukawa term -Fy [a. A, pj, dj, {ai + oq)], in which a is taken as an 
additional fit parameter. The Yukawa term is a function of the densities pj and thicknesses 
dj of the gold, copper, and quartz comprising the Casimir test plates ||T8[| . 

For each value of A investigated, a is then fixed at various values within a range of the fit 
value, and the fit is performed again. The la acceptance region in a is then taken to be all 
the points in the {a, A) parameter space which have a x^/d.o.f. within 1.0 of the minimum 
value. This region, bounded above by the curve in Fig. |1], is roughly one order of magnitude 
more stringent than the Derjaguin limit over the entire range shown, for positive values of 
a [T^. It is very important to note that the measured force plot presumably represents less 



than one percent of the total data set pO[. It is therefore hkely that the hmits attainable 
from a complete re-analysis of this experiment are at least another order of magnitude more 
sensitive than those presented in Fig. |I[ 



C. Other Limits 



Limits on new forces in the range considered here have been obtained in experiments on 
other systems, but they are not sensitive enough to compete with the limits shown in Fig. |l|. 
Bracci, Fiorentini, and Tripiccione |2l| set limits on long-range nucleon-nucleon interactions 



from a re-analysis of a variety of experiments. In addition to macroscopic systems, they 
consider the spectroscopy of hydrogenic and mesic molecules, antiprotonic atoms, and data 
from low energy nucleon-lead and nucleon-nucleon scattering experiments. The limits on the 
coupling strength derived from these elementary systems would have to improve by many 
orders of magnitude in order to appear in Fig. |^. Leeb and Schmiedmayer have derived 
similar constraints from neutron scattering experiments below 10 keV, and from neutron 
optics experiments Studying the difference between scattering lengths measured in 

these two classes of experiments, they extract limits on the size and range of an additional 
Yukawa term in the neutron-nucleus interaction. While more sensitive than those obtained 
in Ref. |2l|, these limits are again many orders of magnitude too weak to appear in Fig. |l]. 



In general, since the electrostatic force between protons is 10^^ times stronger than 
the gravitational force, no experiment on an elementary system which is dominated by 
electromagnetic effects can be expected to provide limits strong enough in this range to 
compete with macroscopic measurements p3| . 



III. THEORETICAL PREDICTIONS 

As indicated in Fig. 0, current experimental limits allow for the existence of forces many 
orders of magnitude stronger than gravity in the range below one millimeter. Predictions 
of new forces — usually involving very light bosons — have arisen in many contexts, from the 
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early notion of Lee and Yang of a field coupled to baryon number to modern ideas 
connected with string theory. For experiments below one millimeter, several well-motivated 
predictions can be addressed. 

Many modern predictions arise in theories which attempt to describe gravitation and 
the other fundamental interactions in a single framework. The mass of the light boson is 
typically given by [Q: 

where Mp ^ 10^^ GeV is the Planck mass (or string scale), M is a mass scale connected with 
the known elementary particles, and is a small integer. For n = 2, the choice of M = 10 
TeV yields m = 10~^ eV, corresponding to a Yukawa force with a range of A = 20/im. 

An interesting class of predictions arises in superstring theories in which supersymmetry 
breaking occurs at low energies (less than about 100 TeV). These theories are of interest as 
a way to address the supersymmetry flavor problem. Superstring models generally contain 
gravitationally-coupled scalars called moduli which are massless at the string scale. Di- 
mopoulos and Giudice consider a "gauge-mediated" model where the moduli acquire mass 
only from supersymmetry breaking, which they suppose to occur at 10-100 TeV P5|. This 



corresponds to Eq. ^ with M = 10-100 TeV and n = 2. The predicted strength and range 
of the modulus forces associated with the down quark, strange quark, and the gluon are 
shown in Fig. |l]. Each prediction covers a region in the {a, A) plane because of uncertainties 
in the model parameters. The predicted strengths are larger than gravity by several orders 
of magnitude, and the ranges are several orders larger than would be expected from Eq. ^. 
Similar scenarios have been discussed previously by Banks, Berkooz, and Steinhardt ||26|| . 

Antoniadis, Dimopoulos, and Dvali have presented a related theory in which supersym- 
metry breaking occurs via Scherk-Schwarz compactification at the weak scale (M ^ 1 TeV) 
p7| . The authors show that theories involving this mechanism are likely to give rise to at 
least one scalar particle, the radius modulus, with a Compton wavelength corresponding 
to the size of the compact dimension. In the model considered, the radius modulus has a 



Compton wavelength in the sub-milhmeter range, and mediates a force with strength 1/3 
that of gravity, as shown in Fig |I[ 

Other hght scalars with stronger couphngs are also possible. The dilaton of string theory 



couples to matter with a strength roughly 80 times that of gravity p8| , p9| . The mass of the 
dilaton depends on unknown physics. Dimopoulos and Giudice suggest a dilaton mass of 
about 10 eV for their scenario, but other authors take the dilaton mass to be an unknown 
parameter The existing limits shown in Fig. Q exclude a mass less than 10~^ eV. 



The axion is a light pseudoscalar boson motivated by the strong CP problem of the 
standard model. Laboratory experiments and astrophysical bounds have left an allowed 
window for the axion mass of 10~^ to 10~^ eV, or A = 20 cm - 0.2 mm [^. However, 
the coupling a of the associated long-range force between unpolarized nucleons has been 
constrained, by recent measurements of the neutron electric dipole moment, to the region 
of the (a. A) plane shown in Fig. |l]. This is probably outside the reach of the experiment 
described in Sec. fV|, but could be a possible goal for the future |^. Moody and Wilczek 
PB| have argued that centimeter-scale experiments with polarized bodies would be more 
promising, and several attempts in this direction have been made P^-^T . 



Still other models are motivated by cosmological observations. The oscillating-G cosmol- 
ogy of Steinhardt and Will involves a massive Brans-Dicke field which is left oscillating at 
high frequency after inflation [^. Yukawa corrections to the Newtonian potential are pre- 
dicted, and useful bounds on a product of the parameters of the model are found when these 
predictions are compared to experiments. The experiment described below would improve 



these bounds by a factor of 600 A recent model due to Beane is based on more general 
assumptions. The author argues that if all experimentally observable phenomena can be 
described in terms of effective field theories, then the observed smallness of the cosmological 



constant implies the existence of new quanta |H0[. Furthermore, if a substantial component 



of the energy density of the universe in the present epoch is due to vacuum energy, then the 
range of the new quanta should be approximately 100 /im, as shown in Fig. |^. R. Sundrum 
develops these ideas further and proposes a specific mechanism that could provide a solution 
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to the cosmological constant problem [^. In this theory, the new quanta are comprised 
of hght gravitational strings. The size of the cosmological constant is set by the scale of 
these strings (as opposed to the standard model scale), below which the behavior of gravi- 
tational forces is expected to be quite different. Consistency with current cosmological and 
laboratory constraints suggests this scale to be between 10 /im and 1 cm. 

Finally, we mention a theory which does not necessarily make recourse to new light 
particles. In a recent discussion, Arkani-Hamed, Dimopoulos, and Dvali propose a theory in 
which gravitational and gauge interactions are united at the weak scale, which is taken to be 
the only fundamental short scale in nature . This model does not rely on supersymmetry 
and is of interest as a solution to the hierarchy problem. The relative weakness of the 
gravitational interaction is a consequence of n > 2 new compact dimensions which are 
large compared to the weak scale, and in which gravitons propagate freely but the standard 
model fields do not. For the smallest number of new dimensions yet to be ruled out by 
experiment {n = 2), the size of the new dimensions is in the millimeter region. Gravitational 
measurements in this regime would be sensitive to potentials dictated by Gauss' law in (4+n) 
dimensions, or gravitational forces varying as r~^. 

From this discussion it is clear that forces of gravitational strength are playing an in- 
creasingly important theoretical role, especially since the advent of string theory. There is 
ample motivation to push to the shortest ranges possible, where macroscopic experiments 
might provide a rare window on physics at the Planck scale. 



IV. SHORT RANGE MEASUREMENT WITH MECHANICAL OSCILLATORS 

As the length scale L of an experiment decreases, the gravitational attraction decreases 
as L^. It therefore becomes difficult at shorter distances to distinguish the Newtonian force 
between the test masses from the background forces. Surface forces are especially trouble- 
some, because they increase rapidly at small separations. The most sensitive experiments 
in the centimeter range to date have used torsion balances. Klimchitskaya, et al., have 
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proposed to use a torsional pendulum inside a rotating steel sphere with a non-concentric 
cavity, and they expect to improve the current limits by roughly an order of magnitude in 



the range above 100 /im [Q. 

Another approach involves the use of small mechanical oscillators for the experimental 
test masses. The experiment of Carugno, Fontana, Onofrio, and Rizzo measured the force 
between a small copper mass and a steel resonator driven at 30 and 1.66 x 10^ Hz, respectively 



4J] . Using a heterodyne detection technique, this experiment was able to set limits on new 
macroscopic forces between, for example, a = 1 x 10^'^ at A = 2 x 10~^ m and a = 8 x 10^ 
at A = 2 X 10~^ m. These limits are roughly two orders of magnitude weaker than those 
shown in Fig. |I|, however, the authors hope to improve their sensitivity by as much as five 
orders of magnitude using fiber optic techniques. 

Our own design is illustrated in Fig. ^ Small mechanical resonators are used for both 
the source mass and force detector. The source mass is driven at the resonant frequency 
of a detector torsional mode, nominally 1 kHz. In the resonant mode of interest, the two 
rectangular sections of the detector counter-rotate. The amplitude of motion is large for the 
smaller section immediately above the source mass, while the larger section is nearly sta- 
tionary. This double torsional arrangement provides a degree of isolation from the oscillator 
mounting, which decreases the mode damping |^ . 

The principal advantage of this high-frequency resonant method is that seismic and 
other vibrational backgrounds can be easily filtered in the kilohertz range. Vibration iso- 
lation is provided by supporting the masses from stacks of brass disks connected by steel 
wires [|^,|^. Other backgrounds include electrostatic forces, which are suppressed with 
a gold-plated shield between the masses, and magnetic forces due to both external fields 
and contamination. Background due to the Casimir effect is comparable to electromagnetic 
effects only for mass separations below 1 fim, roughly 1% of the anticipated minimum gap. 

If these backgrounds can be sufficiently reduced, the dominant source of noise will be 
thermal noise due to dissipation in the detector oscillator. Based on this limitation, we 
calculate the expected sensitivity of the experiment to a new Yukawa-type interaction. 
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The Yukawa force between parallel plates separated by a distance d{t) is given by: 

Frit) = 27iaGpsPdAX^exp {-d{t)/\)[l - exp {-ts/\)][l - exp (-t^/A)] (5) 

where ps and pd are the source and detector mass densities, and ts and td are the thicknesses. 
This would be an exact expression if either plate had area A and the other had infinite area, 
but for the real geometry there are small edge corrections. Neglecting these effects, and 
assuming the geometry in Fig. ^b, the torque exerted on the detector (the smaller rectangle 
only) via a Yukawa interaction with the source is given by: 

Nvit) = 7raGp,pdAdRX^exp{-d{t)/X)[l-exp{-t,/X)][l-exp{-td/X)], (6) 

where Ad is the area of the detector above the source mass, and R is the distance between 
the edge of the detector and the torsion axis. 

The read-out signal from the transducer varies with the amplitude of the applied torque. 
The source mass is driven at the detector resonant frequency ujq, and only the amplitude of 
the signal torque associated with this frequency is effective in driving the detector. From 
Fourier analysis, this amplitude is: 

INyMI = 2TTaGpsPdAdRX''liido/X) exp (-(J/A)[l - exp (-ts/A)][l - exp (-td/A)], (7) 

where Iild^/X) is the modified Bessel function. Here we have used d{t) = d+do sinujQt, where 
do = {draax — dmin)!'^ IS the amplitude of the source mass motion and d = {dmax + dmin)/'^ 
is the average source-detector gap. 

The rms thermal noise torque about the central axis is found from the mechanical Nyquist 
theorem to be: 



AkT (m 



R^\ 



(8) 



\ r \ I 

where m is the mass of the small rectangular section of the detector oscillator, Q is the 
detector quality factor, and r is the measurement integration time. 
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The signal-to-noise ratio is then simply the ratio of equations |^ and ^. Setting the SNR 
equal to unity and solving for a, we obtain the following result for the limit which the 
experiment may set on the Yukawa strength for each value of the range: 

The efficiency factor e, which is of order unity for an optimized geometry, is given by: 

e = 2/i(rfo/A)exp {-d/X)[l - exp (-ts/A)][l - exp (-U/X)] (10) 

The bold curve in Fig. |l|, representing the expected sensitivity of the experiment, is 
calculated using Eqs. |^ and |1^ and the parameters in Table |. For each value of A, the 
maximum separation d^ax is chosen so as to maximize the efficiency e. The maximum 
separation is allowed to vary between the anticipated upper limit of 1 mm and a lower limit 
of 100 /im, below which the backgrounds are expected to dominate. From the limit curve in 
Fig. |1|, we expect the experiment to cover roughly 10 square decades of the {a, A) parameter 
space previously unexplored below 1 mm. 



V. CONCLUSION 

Torsion balance experiments to the present have been sensitive to gravitational-strength 
forces down to millimeter distance scales. While very precise Casimir force measurements 
have recently been carried out and set new limits on effects in the range below 100 fim, 
our calculations reveal that these experiments still allow for forces several million times 
stronger than gravity in this region. Experiments of gravitational sensitivity in this range 
have become all the more imperative with several recent theoretical predictions of new forces, 
motivated especially by advances in string theory. The experiment described in this paper, 
which uses high-frequency mechanical oscillators for the test masses, is expected to reach 
gravitational strength down to approximately 50 /im, and to be sensitive to much of the 
parameter space for these new predictions. 
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FIG. 1. Parameter space for Yukawa-type forces in which the strength relative to gravity (a) 
is plotted versus the range (A). Limit curves from experiments defining the excluded region are 
shown (bold lines) along with the anticipated sensitivity of the proposed experiment, as well as 
specific predictions for new long range forces (fine and dashed lines). 

FIG. 2. a) principal components of the apparatus; b) cross-section of central region showing 
dimensions used in calculation of experimental sensitivity. 
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TABLE I. Parameters used for experimental seusitivily calculaiiou. 

Parameter Value 

Signal Frequency, u/2Tr 1 kHz 

Quality Factor, Q 1 x 10^ 

Minimum Source-Detector Gap, dmin 100/im 

Integration Time, r 1000 s 

Detector Area, 1 cm^ 

Source Thickness, tg 200/im 

Detector Thickness, td 200/im 

Source Density (W), ps 19.3 g/cm^ 

Detector Density (Si), pd 2.3 g/cm^ 

Temperature. T 4K 
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J. C. Long, Experimental Status of Gravitational-Strength Forces in the Sub-Centimeter 
Regime, Figure 1. 
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J. C. Long, Experimental Status of Gravitational-Strength Forces in the Sub-Centimeter 
Regime, Figure 2. 
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